Crossover from coherent to incoherent dynamics in damped quantum systems 
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The destruction of quantum coherence by environmental influences is investigated taking the damped harmonic oscillator and 
the dissipative two-state system as prototypical examples. It is shown that the location of the coherent-incoherent transition 
depends to a large degree on the dynamical quantity under consideration. 
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The destruction of quantum coherence by dissipative 
influences continues to be a problem of central interest 
in atomic physics, condensed matter physics, and chem- 
ical physics. It is relevant to phenomena as diverse as 
wave packet dynamics in atoms |l|, defect tunneling in 
solids , electron transfer in chemical and biological re- 
actions 1^, or quantum computers Q. Important in- 
sight into the effects of an environment on quantum co- 
herence was obtained in the last decade based on stud- 
ies of simple quantum systems in contact with a heat 
bath H,^. Very recently, focusing on the archetypical 
two-state system, two groups |^,p arrived at a result 
which, at first glance, seems to invalidate earlier stud- 
ies. They found that quantum coherence is destroyed 
at significantly smaller damping strengths than thought 
previously, typically differing by a factor 2/3. 

In this Letter, we show that the coherent-incoherent 
transition depends on the particular dynamical quantity 
under consideration (e.g., correlation function, occupa- 
tion probability, etc.). Since different dynamical quanti- 
ties may be associated with different initial preparations 
of the system, quantum coherence may be more or less 
sensitive to dissipation. The resulting critical value of the 
damping strength then changes to a surprisingly large de- 
gree with the respective coherence criterion. This will be 
explicitly demonstrated for the two fundamental dissi- 
pative quantum systems, namely the damped harmonic 
oscillator and the dissipative two-state system. As the 
differences are most pronounced at zero temperature, we 
confine ourselves to this limit. 

We start with the exactly solvable case of a harmonic 
oscillator subject to ohmic damping |Q. The respective 
classical equation of motion for the position q{t) is 



q{t)+7<iit)+^l<lit)=F{t)/M , 



(1) 



where 7 is the usual ohmic damping rate, ujq the fre- 
quency of the bare oscillator, and M the mass of the 
particle. The response to the external force F{t) is de- 
scribed by 
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dt'xosc{t-t')F{t') , 



(2) 



where {q(t)) denotes the expectation value, and Xosc{t) 
is the linear response function, which for the harmonic 
oscillator coincides with the classical response function. 
One dynamical quantity of interest is the equilibrium 
correlation function 



asc(i) = Re(g(i)g(0)), 



(3) 



where the brackets denote a ground-state average. Con- 
sider now the associated spectral function 0] 

S'osc(w) = Cosc{t^)/\t^\ , 

which at T = is related to the absorptive part of the 
dynamical susceptibility by (we put Ti = 1) 

S'osc(w) = ImXosc('^)/w ■ 

It is an even function of lo, and the second form follows 
from the fluctuation-dissipation theorem. From Eq. (|^), 
the dynamical susceptibility 



( \ = ^0 

and hence the spectral function 
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(cJq — uj'^Y + 7^0;' 
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(4) 



(5) 



Another problem of interest is the relaxation of the 
expectation value {q{t)) starting from a nonequilibrium 
initial state. Applying the force F{t) = Afu;Qqo8(— i), 
the initial condition (q(0)) — qo is prepared [cf. Eq. (||)], 
and the relevant dynamical quantity is 



Poscit) = {q{t))/qo 



(6) 



With the Fourier transform (Q) of the response function, 
one obtains from Eq. (0) 



Posc{t) — cos{ilt — (j)) exp(— 7^/2)/ cos(0) 



(7) 



where fi = -^/wq — 7^/4 and (/> = arctan(7/2J7). 

Let us now discuss the coherent-incoherent transition 
as a function of the dimensionless damping strength 

a — 7/2(jJo • 

We first consider the appropriate coherence criterion 
based on the spectral function S'osc(w). For weak damp- 
ing, a < oic, the function Sosc{^) exhibits two inelastic 
peaks at finite frequency oj = ±Wm(a). At the critical 
damping strength Uc, one has LOm{cic) = 0, and the two 
peaks merge into a single quasi-elastic peak centered at 
a; = 0. This quasi-elastic peak then persists for a > ac- 
The value of ac can be determined by inspecting the 
sign of the curvature of S'osc(^) at zero frequency. From 
Eq. (pi), we have the small-cj expansion 



Sosci^) = 2axl [1 + (2 - Aa')xl oj^ + 0(c^")] 



(8) 



where Xo = l/i^o is the static susceptibility. Thus the 
curvature of Soac{^) is positive (implying coherence) for 
a < 1/ V2, but changes sign at the critical value. Regard- 
ing the equilibrium correlations Coscit), quantum coher- 
ence is suppressed for a > ac = l/v2. 

A different coherence criterion can be developed based 
on the quantity Posc(t). For weak damping, one finds 
damped oscillations such that Posc(i) changes sign occa- 
sionally. As the damping strength a is increased, Poscit) 
shows a transition from damped oscillatory to a purely 
incoherent behavior. The coherent-incoherent transition 
is reached at the critical value a — a*. For a > a*, 
one has Posc{t) > for all times t. This coherence cri- 
terion based on a noncquilibriunr initial preparation has 
been used in most previous work p, |q|. For the damped 
harmonic oscillator, we see from Eq. (|7|) that Posc(i) be- 
comes overdamped for a > 1, while damped oscillations 
persist for a < 1. Regarding the nonequilibrium quan- 
tity Posc{t), one has destruction of quantum coherence 
for a > a* — 1. The ratio between the critical damping 
strengths for the two dynamical quantities Cosc{t) and 
Posc{t) is then given as 



ac/a* = 1/^2 . 



(9) 



We now turn to the case of a symmetric two-state sys- 
tem. Coupling to an ohmic heat bath is described in 
terms of the spin-boson Hamiltonian p,p| 



H = -(A/2)a, 



E 



2m,- 



CiX^C;, 



The two eigenstates of az with eigenvalues ±1 are coupled 
by the transfer matrix element A representing the tunnel 
splitting of the free system. The ohmic bath of harmonic 
oscillators is fully characterized by the spectral density 



^M-J^ 



where a is a dimensionless damping strength, and lUc is 
a high-frequency cutoff. 

An important aspect of this model is its correspon- 
dence with the Kondo effect. By writing the parti- 
tion function in the Coulomb gas representation, a firm 
equivalence with the anisotropic Kondo model can be es- 
tablished M. This correspondence has been exploited 
by Costi and Kieffer 0, who applied a dynamical ver- 
sion of Wilson's numerical renormalization group to the 
anisotropic Kondo model. Furthermore, Lesage, Saleur, 
and Skorik H have employed integrability and form- 
factor techniques to study this model as well. Both 
groups have calculated the ground-state spin-spin cor- 
relation function of the Kondo model equivalent to the 
T = equilibrium two-state correlation function 

Cit)^Re{a4t)a40)) . 

This quantity is the direct analog of the harmonic oscil- 
lator correlation function Coscit) defined in Eq. (H). 

Another useful quantity, particularly in the context of 
macroscopic quantum coherence [|||§], is the occupation 
probability 

P{t) - {<Jz{t)) , 

corresponding to the harmonic oscillator quantity Pose (t) 
defined in Eq. (g). In contrast to C(i), the function 
P{t) is subject to the nonequilibrium initial preparation 
az(t = 0) = -\-l. This preparation of the initial state may 
be realized by applying a large external bias. Thereby the 
spin is held fixed in the state az — +1 with equilibrated 
environment. At time zero, the constraint is released, 
and the dynamics starts out from P{0) — I with this 
factorized system-environment initial state. 

Let us now discuss the coherent-incoherent transition 
for the dissipative two-state system, starting with the 
equilibrium correlations. Similar to Eq. (||), the function 
S{ll}) = C{u})/\u!\ has the low- frequency expansion 



S{uj) = 2^X0 [1 + «(«)Xo ^^ + 0{uj^)] 



(10) 



where Xo is the static susceptibility and K(a) is a dimen- 
sionless parameter. The zero-frequency limit of Eq. (nfl) 
is the generalized Shiba relation for the spin-boson prob- 
lem |12]. As a result of this relation, spin-spin correla- 
tions decay asymptotically as C{t) = —2axo/t'^- In Refs. 
[n,p|, the function K{a) has been computed and was found 
to change sign at the critical value ac = 1/3. 

Next we analyze the transition from oscillations to in- 
coherent relaxation in the quantity P{t). Within the 
widely- used noninteracting-blip approximation (NIBA), 
the critical value is a* — 1/2 M. This result can easily 
be seen by switching to the Laplace transform P{X) and 
defining a self-energy S(A), 



2 ^—^ miU}; 



■5{uj — LUi) — 2TraLU exp{—uj/LUc) , 



P(A) = 1/[A + S(A)] 



(11) 



NIBA gives for the self-energy p,p 



S(A) = A^ cos(7ra) / dr 



exp(— Ar) 



2a 



/O {^cT) 

= A,(A,/A)i-2a^ 

with the effective frequency scale 

Ae = A[cos(7ra)r(l - 2a)]i/2(i-")(A/w,)"/(i-") 



(12) 



In the limit 



1/2, the frequency Ae approaches 



7rA^/2cJc- For a > 0, the function P{X) has a branch 
point at A = 0. The complex A-plane is cut along the 
negative real axis, and in the cut plane P(A) is single- 
valued. The cut leads to an incoherent contribution to 
P{t). Moreover, for a < 1/2, the integrand has a conju- 
gate pair of poles in the cut A-plane describing damped 
oscillations. 

To study the coherent-incoherent transition, wc con- 
sider the pole condition from Eq. (pT|), 



(A/A< 



\2-2a 



-1 



(13) 



and put 



a = l/2-e, |e|<l. (14) 

For e > (a < 1/2), insertion of the ansatz 

A/Ae = -1 + e-u ± iev + 0{e'^) 

into Eq. (O) yields u = and v — 2tt. Hence, in the 
NIBA the damping rate F and the oscillation frequency 
il are up to corrections of order e^ 



F/Ae = 1 , 0/Ae =: 27re 



(15) 



This results in damped oscillations with frequency Q. In 
contrast, for e < (a > 1/2), the poles are not in the cut 
plane and therefore give no contribution to P{t). In that 
case, P{t) is completely given by the incoherent branch- 
cut contribution. Thus the critical damping strength is 
indeed a* = 1/2. 

Remarkably, exactly at the special value a = 1/2, 
NIBA becomes exact [gj, while it is only an approxima- 
tion for a 7^ 1/2. One serious deficiency comes from 
the branch-cut contribution which would imply the exis- 
tence of an algebraic long-time tail. From Eq. (f2h, for 
a = 1/2 - e and A ^ 0, one finds P(A) ~ X^". The 
asymptotic branch-cut contribution is therefore 



Pit) = -2e/(Aei)i- 



-2e 



(16) 



This term decaying slower than 1/i^ contradicts the 
fluctuation-dissipation theorem. Hence, the long-time 
tail (^6|) is an unphysical artefact of NIBA g] . Such a fail- 
ure raises the question whether the NIBA value a* — 1/2 
for the coherent-incoherent transition remains correct. 



To investigate the coherent-incoherent transition be- 
yond NIBA, we now systematically expand around the 
exactly solvable case a = 1/2. In diagrammatic terms, 
the exact self-energy S(A) is the sum over all irreducible 
arrangements of "blips" and "sojourns" . A blip (sojourn) 
refers to the time spent in an off-diagonal (diagonal) state 
of the reduced density matrix MM- For instance, the 
NIBA expression (021) is just given by the single-blip con- 
tribution, thereby effectively disregarding all inter-blip 
interactions in P(A). The fact that NIBA becomes ex- 
act for a = 1/2 is explained simply by the concept of 
collapsed blips P]. In view of Eq. (|l4h, there is a factor 
cos(7ra) = 7re in Eq. (|l2|). This 0{e) factor must be can- 
celled by a 1/e short-time contribution of the r-integral 
over the length of the blip such that a finite result can 
arise. Therefore, for a -^ 1/2, only blips of effectively 
vanishing length, that is "collapsed" blips, contribute. 
Since interactions among different collapsed blips vanish, 
NIBA becomes exact [||. 

For nonzero e, the blip length is finite, and one has 
to take into account all sequences of collapsed sojourns 
within an extended blip of length r. Their presence cru- 
cially modifies the self-energy in the limit A — > 0. Since 
collapsed sojourns do not interact with blips or other 
sojourns, a grand-canonical gas of collapsed sojourns 
merely gives a factor exp(— Aer/2), and the self-energy 
13) is changed into 



/>OC 

E(A) = AVe / dr 
Jo 



CXp[-(A + Ae/2)T] 



Ae(A/Ae + l/2)-2^ . 



(17) 



With the regularized self-energy (|l7|), the pole condi- 
tion now reads 

A(l/2 + A/Ae)'^--Ae. 

This yields the exact decay rate and oscillation frequency 
up to corrections of order e^ , 

F/Ae = 1 - 2e In 2 , Q/Ae = 2iTe . 



Compared to the NIBA result (|15|), the rate acquires a 
correction in order e. However, the oscillation frequency 
remains completely unchanged. Therefore, the NIBA 
value a* — 1/2 for the location of the coherent-incoherent 
transition turns out to be the exact result. This implies 
the critical ratio 



ac/oi*c = 2/3 , 



for 



which is slightly smaller than the respective ratio 
the damped harmonic oscillator. 

The exact self-energy (O) shifts the branch point of 
P{\) and thereby removes the spurious algebraic long- 
time tails (|l^. For e > 0, the leading branch-cut contri- 
bution at long times is now given by 



P{t) = -2e 



exp(-AeV2) 



while for e < 0, we obtain 



Pit) - 8|6| 



exp(-Aet/2) 



Thus, the unphysical algebraic long-time tails are sup- 
pressed by an exponential decay factor. It is straightfor- 
ward to see that for e > (a < 1/2), the full cut con- 
tribution is negative, while it becomes positive for e < 
{a > 1/2). Indeed, for a > 1/2, the function P{t) is 
positive and monotonically decaying, i.e., the dynamics 
is fully incoherent. In marked contrast to the NIB A re- 
sult dlq), the power of the algebraic decay factor does 
not depend on the sign of e. 




FIG. 1. Monte Carlo data for P{t) at a = 0.4 and T = 0. 
Circles are data points for A/uuc = 1/6, the solid curve con- 
necting them is a guide for the eye only. Statistical errors are 
well below 5%. The dashed curve is the NIB A result. 

The remarkable success of NIBA in predicting the cor- 
rect value of a* provokes questions about the general 
quality of such a simple approximation for the full range 
of a. A convenient tool to investigate this issue is the 
real-time quantum Monte Carlo simulation method pi . 
This technique permits a numerically exact calculation 
of both dynamical quantities C{t) and P{t) by stochastic 
evaluation of the respective real-time path-integral rep- 
resentations. The dynamical sign problem arising from 
the interference between different real-time paths can 
be largely circumvented by a partial summation scheme, 
and stable simulations can be carried out for rather long 
times. The simulation code permits a computation of 
P{t) directly at zero temperature, where the case a = 1/2 
serves as a convenient benchmark which is indeed passed 
accurately [ p| . 

Numerical results for P{t) at a = 0.4 are shown in 
Fig. n^, and data for a > 1/2 can be found in Ref. pi[ . 
These data represent universal scaling curves in the sense 



that they can be obtained from different ujc and/or A by 
employing the effective frequency scale Ag. Apparently, 
on short-to-intermediate time scales, NIBA yields a very 
accurate prediction for P{t). However, the prediction 
C(t) = P(t) as well as the long-time tails dlq) are clear 
deficiencies of this approximation |11|. 

In conclusion, we have demonstrated that the initial 
preparation of a dissipative quantum system leads to 
drastic changes regarding the transition from coherence 
to incoherence as the damping strength is increased. 
Both the damped harmonic oscillator and the dissipa- 
tive two-state system show that the two coherence cri- 
teria previously employed in the literature lead to a 
factor l/-\/2 and 2/3 difference in the critical damping 
strengths, respectively. Any investigation of the environ- 
mental destruction of quantum coherence thus necessi- 
tates a clear specification of the physical quantities under 
study. 
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